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Abstract 

We present the way the Lorentz invariant canonical partition function for 
Matrix Theory as a light-cone formulation of M-theory can be computed. We 
explicitly show how when the eleventh dimension is decompactified, the M = 1 
eleven dimensional SUGRA partition function appears. We also provide a high 
temperature expansion which captures some structure of the canonical partition 
function when interactions amongst D-particles are on. The connection with the 
semi-classical computations thermalizing the open superstrings attached to a D- 
particle is also clarified through a Born-Oppenheimer approximation. Some ideas 
about how Matrix Theory would describe the complementary degrees of freedom 
of the massless content of eleven dimensional SUGRA are discussed. Comments 
about possible connections to black hole physics are also made. 
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1 Introduction 



It has been intense the attention that the proposal known as Matrix Theory |] 
has attracted from the community of physicists working on the field String Theory 
has become after the so called second string revolution. It seems that after a first 
boost of frenetic activity few physical quantities have been computed so as to 
test to what extent we can at least recover what we knew before from D-brane 
dynamics. At present, it could be considered pretentious trying to get exact non 
perturbative results from Matrix Theory. In Theoretical Physics, the way out 
of this kind of situations is to wait for a more powerful formulation of what we 
have at hand. Today's form of Matrix Theory has some drawbacks. First of all 
the original formulation was one in the IMF. This means that one should take 
a subtle N — > oo limit to recover Lorentz covariance. A reformulation of the 
conjecture by using discrete light-cone quantization (DLCQ from now on) seems 
to provide a meaning for the computations at finite N. In DLCQ, it is the light- 
like coordinate x~ that gets compactified. In this form, the conjecture would 
establish a correspondence between finite iV Matrix Theory and the DLCQ of M- 
theory. Support for this correspondence has been given in ||. Another related 
point is that the present form of Matrix Theory is not a background independent 
one. A list of Matrix Models for given backgrounds is in need and, up to now, 
things do not go very fast to this respect. Something as modest as a formulation 
of Matrix Theory for toroidal backgrounds is still lacking. Consequently, Matrix 
Theory might be nothing more than another step towards getting a more complete 
description of what in the past we called String Theory. It is then not clear how 
far we will be able to go. 

It seems however that there are some ingredients that could show up as ad- 
vantages. After all we have a Hamiltonian and perhaps this could be used to 
compute physical quantities that from their own definition are adapted to that 
formulation. This would be the case of the canonical partition function, Z(f3). 
This work will be devoted to computing this quantity, although the task will not 
be exactly accomplished. The interest of the work will be that of getting a better 
knowledge of what Matrix Theory Thermodynamics might be. From the very 
beginning it seems clear that with T — as the thermodynamical temperature 
we use to define the canonical ensemble, the canonical partition function appears 
to have only sense in the DLCQ picture because the IMF formulation of the con- 
jecture is actually throwing away energy configurations it seems one should take 
into account at finite temperature. We will focus on the case in which all spatial 
dimensions are open with the exception of x~ . One might argue that owing to the 
fact, also known in string perturbation theory, that compact spaces increase the 
number of degrees of freedom even to the extent of making difficult a formulation 
of a Matrix Model for T n for n > 6, one should compute Z (/3) for a T 9 space 
and get from here Z (j3) for R 9 . This prevention assumes that this tremendous 
amount of degrees of freedom do not fully disappear after decompactification. 
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Such an attitude is based on the extended misconception that a microcanonical 
description of perturbative string theory gives a decompactification limit which 
is different from that of the microcanonical picture gotten directly from Z (/?) for 
an already open box (cf. ||). Then there is no a priori reason to think that a 
formulation of Matrix Theory on a T 9 space, if there is such a thing, would not 
give what we have now as its decompactification limit. In fact, by computing 
the free energy of the D-strings, by the way, we will check that the extra degrees 
of freedom from the extra compact dimension decouple when its length goes to 
infinity. 

For the time being, in section 2, we will expose the technique to compute 
Z (P) for a theory which is formulated in the light-cone frame. The main point 
is that if T = 1/(3 can be treated as the temperature one could measure with a 
thermometer, one has to admit that it has to transform as the energy, i.e. p°. 
In other words Z ((3) = Tre _/3p is a Lorentz invariant quantity we are going to 
compute knowing only a description in a light-cone frame. We will present this 
computation in the DLCQ form of the Matrix Theory conjecture. In section 3 
we will study the limiting case in which one has free (classical geometry) con- 
figurations and what happens when off-diagonal terms (open stretching strings) 
are taken into account. In section 4, we will study a classical system of matri- 
ces or what more properly could be defined as the classical statistical mechanics 
associated with Matrix Theory and we will relate it to Matrix Theory at finite 
temperature by computing its quantum corrections. In section 5 we will pertur- 
batively compute the partition function by a finite temperature extension of the 
computation in ||. We will finally get out the conclusions in section 5. 

2 The canonical free energy from a light-cone 
frame description 

The canonical partition function for the simplest relativistic system of a gas of 
particles of mass m is by definition (see M) 



This is an invariant quantity whenever we admit that the temperature T trans- 
forms as the energy, i.e. the zero component of the momentum. This is in physical 
terms a consequence of the interpretation of the temperature as measuring the 
average kinetic energy of the system. If the system is known in the light-cone 



Z(j3) = Tre-W. 





Z(0) 
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Let us perform the trace by first taking a very illustrative path. One can 
compute a single particle partition function by performing the trace in (^) using 
a basis \p + > <8> \p T > of single particle momentum eigenstates. If the space 
volume goes to infinity, we have 

L f+°° 0p+ m2 P L dr f > _^fz 

-I ip * e -^ e -^—^j^ pe -^. (3) 

The change of variables s = and an elementary integration take us to 

T / ds s 2 e 2s e 2 . (4) 

(2tt)2 Jo 

Where d is the space-time dimension, d — d T + 2, and L — ► +oo. The vari- 
able s = is the Schwinger proper time, and we have the well known proper 
time representation of the Helmholtz free energy. Eq. (Q) can be recognized as 
—/3Fi (/?) = Z\ (j3) for a particle of mass m and also as the same magnitude per 
degree of freedom for a quantum field with particle excitations of mass m. In 
both cases one has Maxwell-Boltzmann statistics. 

To get Z(f3), i.e. the canonical partition function, a further exponentiation 



and another sum are needed (see |10|). The relation, for particles obeying Bose 



statistics, is 

m^g^M. (5) 

The exponentiation gives the multi-particle system and the sum over r gives 
bosonic quantum statistics versus Maxwell-Boltzmann which corresponds in fact 
to the r = 1 term in the sum. So one has two aspects which, in principle, are 
different. One can have a multi-particle description but with Maxwell-Boltzmann 
statistics instead of Bose or Fermi quantum statistics. This subtlety shows up 
clearly as soon as one deals with the computation of the canonical partition func- 
tion for string perturbation theory in the light-cone gauge. It is worth noticing 
here that the single-particle partition function depends linearly on the open (infi- 
nite) volume and then the multi-particle partition function depends exponentially 
on the volume. 

Eq. (H) is the result of computing Z {(3) as the trace of e _/3p using a basis, 
in the Fock space, of multi-particle eigenstates of the energy p° assuming Bose 
statistics. The rigorous way of getting (|5]) is assuming that our spectrum is 
actually discrete. Physically this can be achieved by enclosing the system in a 
finite volume box. One can then write 

Z{(3) = J{{l±e-^f 1 . (6) 

k 
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Here is the energy, that is usually understood as a dispersion relation and the 
signs distinguish bosons from fermions. Taking the logarithm of both sides, one 
gets the Helmholtz free energy as an infinite sum. It is from here that one gets 
the sum over r by expanding In (1 ± x). Physically it is not necessary to assume 
that the compactified theory is known if one makes the hypothesis that, after 
taking the infinite volume limit that convert the sum over discrete momenta into 
an integral, the remnant degrees of freedom are those of the formulation of the 
system for an infinite box which is supposed to be what is known. 

In quantum field theory one can check that In Z (/?) can also be obtained 
by a Euclidean path integral computation of the covariantly defined theory on 
W 1 ^ 1 x Sp. After clarifying some questions related to modular invariance, the same 
correspondence will work for perturbative string theory, but a light-cone gauge 
computation based on the transverse Hamiltonian does not produce a multi-string 
partition function on its own. Let us explain what we mean. The starting point 
is again Eq. (p]) which can be rewritten as 

ZW = ± rVe^ T (-AJ (7) 
2n Jo \2p + a J 

Here, zt is defined as a transverse partition function given by 

z T (A) = Ti Ht e- XHT (8) 

with Ht the light-cone gauge Hamiltonian associated with the light-cone gauge 
(transverse) action that for the bosonic part reads 

S l h L Bic = J dadrd a X*d a X* i = 1,...,8 (9) 

A is a sort (it is a dimensionless quantity) of inverse transverse temperature. 
With this action (and is fermionic part), Zt can be computed by performing 
a Euclidean path integral with periodic (anti-periodic) boundary conditions in 
the transverse temperature A for the coordinates (32-spinor S) as fields in two 
dimensions. The result one finally obtains for a string one-loop computation is 
again a Maxwell-Boltzmann Helmholtz free energy proportional to the infinite 
volume in nine dimensions; i.e. a single string partition function 2| trmg {(3) [|10|j . 

The supersymmetrical version of Eq. @ gives the free energy for the multiple- 
string gas. In short, it can be written 

\*Z{fS) = £ (10 ) 

r=odd>0 r 

Where one assumes that the single partition function includes the sum over 
Maxwell-Boltzmann contributions coming from the bosonic and fermionic degrees 
of freedom. It is worth noticing again that the single-string partition function de- 
pends linearly on the infinite volume while the multiple-string partition function 
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exponentially increases its free energy with the volume. In terms of the canonical 
entropy, it is clear that the Helmholtz free energy has to depend linearly on the 
volume as the entropy does. 

Let us finally tackle our goal which is to get the form of the partition function 
for the Matrix Model. One first notes that in the DLCQ, p + is discrete because 
x~ is a circle of radius R. From an eleven dimensional point of view, what one 
would first write down for the partition function of the model would be 



+00 

? 



^matrix = & - /3 m rj^, £ - /3RH {N) ^ -Q 

jV=l 

where H = RH^ is the Hamiltonian for the SYM with U(N) symmetry given, 
with the gauge A = 0, by 



h = rtt -n 



1 



r2 



2 lQn 2 a 



Y\Y 3 



2 



2ira' 



3/2 



Trypan (12) 



Physically Eq. ( |TT]) should be interpreted as the classical partition function for a 
single eleven dimensional object. The question is that we have to determine what 
is meant by the trace after the sum. After knowing the features Matrix Theory 
have revealed since the conjecture was formulated, it appears reasonable to us 
to study particular configurations in order to try to identify the single object 
one could use to make statistical mechanics the way it has been presented here 
for particles, fields and perturbative strings. At first sight one knows that in 
terms of ten dimensional physics one has systems of N free DO-branes and also 
bound states at threshold formed by iV DO-branes. Let us then study the simplest 
picture which is that of the free (classical geometry) configurations in the Matrix 
Model. 



3 Matrix Theory free configurations 

A simple analysis of the Hamiltonian we have written down shows that in Matrix 
Theory there is no well defined meaning for the concept of position interpreted as 
a way of defining a standard geometry for the configuration space of the blocks 
we play with which are DO-branes. The reason is well known and results from 
the natural thickness of Dirichlet branes given by the fact that they have open 
strings on them and actually at ultrashort distances stretching strings can glue 
together two DO-branes, for example. These stretching strings actually introduce 
some kind of fuzziness that makes a classical geometrical interpretation for the 
configuration space of the branes an approximate concept. If the N x N matrices 
we use to write H(n) are diagonal, then a great simplification results because 
commutators get null and the Hamiltonian reduces to the simple form 
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Hfc = Tr (In?) , (13) 

with each being a diagonal N x N matrix. It is clear because of the form 
of the Hamiltonian that the classical (geometry) configurations correspond to 
a free theory. The states we want to perform the trace over in Eq. (|TT|) are 
[(®i=ib l >) ® \t >} n with the subscript N indicating that the states are actually 
arranged in an iV component vector for a given N and | r > carries the Pauli spin 
information. With them, one can make the following guess 



+ 00 . +OO I T/JV / 0_ \ 2 N 

^£ e -^ e -^ = £^-^(g e -^. (14) 

Here the N\ factor comes from the permutation group which is the discrete sym- 
metry that survives from U(N) after truncating the matrices to the diagonal 
ones. Vg is the spatial nine- dimensional volume. 

Now the next step would be that of substituting this single partition function 
into Eq. flT0|) and getting Z (/?) after an exponentiation. However there is some- 
thing remarkable in Eq . fll4]) . It has been constructed as a single particle function 
an as such should be identified with —(5F\ {(3) and then one would expect it to 
depend linearly on the volume and this is not the case. In fact one can write 



Z l { cc = exp 



256 




-1, (15) 



256 7^9(^1 e~ p ™. (16) 



i.e., this partition function depends exponentially on the nine-dimensional vol- 
ume. This is simply the result of the fact that the Matrix Model contains second 
quantization in itself in the sense that the first quantization of an object in eleven 
dimensions gives the description of ten dimensional multi-objects which are sets 
of free DO-branes. Looking at it this way, Eq. ( [[5f) would define the single object 
as one with single partition function 

V9 
(2tt) 1 

The second quantization of this single object would produce composite systems 
of free DO-branes. In other words, a given total momentum would result from the 
sum of the individual momenta of each free DO-brane. Then the special bound 
states at threshold, special because they are bound states of n DO-branes but 
with zero relative energy, are missing and we are not taking into account that the 
given total momentum can also be shared by several bound states at threshold. 
The consequence is that we have to modify our single object to include these 
states and that is easy. Our single object will be a sum over partons of any 
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positive RR charge. This is similar to the image of the single fundamental string 
in the analog model as a collection of an infinite number of fields with masses 
running from zero to infinity. At last, the single object partition function will be 



+00 



Z\ cc ((3) = 256 



Vo 



g V 2R 



k=1 (2nf \ R/3 , 
With this, we can now apply the recipes in Eq. ([IT]) to get 



{(3) = exp 



1 



_r=odd>0 



(17) 



18) 



The relationship between Z bee (/3) and our starting point in Eq. ( |i~4| ) is very 
clear. Eq. (14) coincides with the r = k = 1 term in Z iicc {(3) after expanding the 



exponential. In other words it is the multi-object Maxwell-Boltzmann version of 
our system without counting the DO-brane bound states at threshold. 

We see in ( |TT| ) that our single object is a superposition of bound-states com- 
posed by k D-particles. The dynamical degrees of freedom of each parton does 
not contain the dynamics of its components, i.e. we do not have Vg or any per- 
mutation property of the k D-particles as in a standard gas with the same number 
of components; it only contains the dynamics of the center of mass. This feature 
is interesting because it does not prevent the interpretation of a given fc-parton 
as a Schwarszchild black hole (see ||). We will explain it in more detail in section 
5. 

It is necessary to check the nature of our single object by identifying into what 
the single partition function turns after taking the R — > +00 limit that would 
open the eleventh dimension up. It is an easy task to perform such limit because 
amounts to converting the infinite sum over positive integers k into an integral. 
One has 



2tt Jo (2tt) 9 \ (3 j 

2567r- 11 / 2 r(^) F 10 /r 10 , (19) 

with V10 the space volume in eleven dimensions. This, through Eq. (|i~8D, gives 
— (3F (f3) for the massless field content of M = 1 SUGRA in eleven dimensions. 
Our single object becomes an eleven dimensional supergraviton when opening 
up the light-like dimension as, may be, one could have expected from the very 
beginning. At finite temperature we have to sum over the complete tower of 
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longitudinal momentum modes, that is, we have to take into account the contri- 
bution of all the finite-N Matrix Models, this allows us to avoid the problem of 
relating the finite-N model to the eleven dimensional supergravity [|12|]. 

The question now is whether we can use all this analysis to treat interactions 
among the DO-branes. At first sight one would say that all the treatment for 
particles, fields, strings and the Matrix Model itself is only accurate for free 
single objects. We know this is not true for fields and strings. The proper 
time formalism is also accurate for including interactions perturbatively through 
corrections 5m? to m 2 . One only has to change the mass squared by the mass 
squared plus the loop corrections. This can be done for strings in the analog model 
in which the free energy is the sum over the field content of the string (cf. |[21|| ). 
We will actually see in section 5, that this can also be done for the Matrix Model 
in the approach in which a perturbative expansion around a fixed background is 
done. In some sense, what the light-cone description of M-theory should provide 
for the massless content of Af = 1 SUGRA is analogous to the massive tower of 
ten dimensional Planck masses that promote the ten dimensional supergravities 
to the vibrational modes of the full fundamental superstring theories. The Matrix 
Model interactions through stretching strings are the sources for these massive 
companions, although we do not see any a priori reason to believe that the relevant 
degrees of freedom could be described as some kind of massive quantum fields 
beyond this fixed background perturbation expansion approach. 

4 Statistical mechanics of U(N) matrices 

We have already stressed as one of the advantages of the present formulation 
of Matrix Theory that it is based upon a relatively simple Hamiltonian. This 
provides us with a different way to do our computations without any need to do a 
hard path integral nor to restrict ourselves to a particular perturbation series that 
strongly depends on the background -the vacuum- we choose to expand around. 
At least, as we shall see, it is possible to find global, non perturbative features of 
the partition function that the field perturbation theory cannot capture. 

The Hamiltonian formulation fundamentally differs from the Lagrangian one 
in the mathematical objects that contain the complexity of the system. In the 
Lagrangian formulation, the function that defines the system is the action, cal- 
culated over the trajectories or paths. This necessarily leads to a path integral if 
we want to calculate any physical quantity. On the other hand, in the Hamilto- 
nian formulation, the system is defined by a hermitian operator acting over wave 
functions. Physical magnitudes are calculated through summations or integrals 
over real or Grassmann variables. The complexity is thus transported into the 
search for appropriate wave functions as well as their matrix elements with the 
operators that appear in the Hamiltonian. 

Let us rewrite the transverse Matrix Theory Hamiltonian with every h and 
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a' in it 




16a' 3 ft 2 L 



1 



Y\ Y j - 



27r(a'h) 3 / 2 



i 



iry> [Yj,6] 



) 



(20) 



with the constraint 



TT = 




(21) 



Again, we have chosen the gauge Aq = 0. Besides, one has to impose another 
condition that arises from the Gauss law. If we take the Arnowitt-Fickler gauge 
in the 9 + 1 theory, the condition after dimensional reduction to + 1 is = 0. 
The only consequence is that now, the bosonic index i takes values from 1 to 8. 
This gauge election does not need ghost fields. This is general for all the 'gluons' 
except for that related to the identity. This one decouples from all the rest and 
does not even appear in the Gauss constraints because trivially commutes with 
any other matrix 

More comments on the particular form of this Hamiltonian are also in need. 
Four different operators appear: the fermionic and bosonic positions 9 and Y 
and their conjugate momenta IT and tt. The commutation relations among them 
are various. Firstly, all of them are N x N matrices and therefore they obey the 
U(N) algebra. Secondly, they are either bosonic or fermionic creation operators 
of one- dimensional world-line fields or, equivalently, operators with either real 
or Grassmann numbers as eigenvalues from the ten dimensional point of view. 
So they correspondingly commute or anti-commute according to their nature. 
Finally, each coordinate and its conjugate momentum have canonical Heisenberg 
commutation rules too. 

Other noticeable characteristic of the Hamiltonian is the ubiquitous appear- 
ance of h in every term. The physical origin of this fact is that, from the eleven- 
dimensional point of view, the BPS masses of the DO-branes are Kaluza-Klein 
energies and, therefore, purely quantum mechanical. Besides, let us remember 
that the Planck length in any dimension has also got a purely quantum me- 
chanical nature. As it stands, it is impossible to say what is the limit of the 
Hamiltonian as % — > 0, or even if there is any. 

The first step we shall take is to expand the matrix operators in terms of the 
generators of the U(N) algebra in order to avoid complications related to their 
matrix nature. This yields 

2 In the special case of a + 1 dimensionally reduced theory in the hamiltonian formalism, 
the Gauss constraints can be completely written in terms of the structure constants of the 
gauge group, that is §I| G p = F^YflUj - iFjfOW. Let us note that in the free case in which 
the gauge symmetry is broken to U(1) N these conditions do not impose any extra constraint to 
the election Aq = and so we keep nine physical scalar fields for each D-particle. This justifies 
the calculation in section 3. 
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D D AT) 
ff — _n^TT^ — F^FTeyiyiyjyj ^ ayi jpcv Q /QO^ 

H ~ 2h VhVh lfora**' p 2n(ha'f 2 7 ^ ^ 

where F° v are the structure constants of the gauge group. To make the trace 
over the exponential of this Hamiltonian we have to choose a base of wave func- 
tions. The most appropriate is formed with eigenvectors of the 16iV 2 fermionic 
coordinates and the 8iV 2 bosonic momenta. We could notate them by |n^0" > 
where the indices mean that the state is defined with the (16 + 8)iV 2 eigenvalues. 
This quantum numbers do not always completely determine the system. At least 
we know that when two or more of the coordinates coincide there appear bound 
states at threshold that cannot be described if we do not add another number. 
We shall ignore this fact since we shall not give precise numerical results in this 
section and this degeneracy would only alter the value of certain coefficients. The 
complete calculation we want to perform is 



z = f e-m f (n dn*) f n do 1 ) < n^|e-^KW) \ U fy > (23) 

N=l J \k,X J J 

We shall first carry out the integral and then sum over the index N. Let us 
begin with the integral over the fermions 

I f = (W) 4 ^ 2 - 1 ) I d9e- eaM « 9 » (24) 

with 

{e a ,6 b } = 6 ab VKa! (25) 

and where the indices shown are labels related to the gauge group as well as 
to the sixteen different fermionic degrees of freedom. The constant before the 
integral comes from the relation between the coordinates and their momenta. Its 
function is to make the integral dimensionless. 

M = iph- 1 a'- 2 R^-Y^F^ li (26) 

It is easy to expand the exponential in a Grassmann series. The matrix M is 
16(iV 2 - 1) x 16(iV 2 - 1) so that the 8(iV 2 - l) th term is the first to contribute. 
Here again, the fermions related to the identity are not taken into account be- 
cause they do not appear in the Hamiltonian. As regards to this first term, 
we can proceed just ignoring all the anti-commutators because any term coming 
from them and depending on h cannot have all the variables we need for the 
integral not to be null. Therefore, it is possible to order the operators separating 
the coordinates from the momenta, and to convert it into a double integral of 
Grassmann numbers, not operators. The result is known to be 
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JdB^ E <e\\{{-e l m l] e 3 )\e>= 

combinations 




The next term does not contribute because there is a repeated coordinate. 
There are some terms with 8(iV 2 — 1) + 4 operators that do contribute. Let us 
show how with an example 

9 1 m 1 29 2 9 2 m2i9 1 [■■■]= m 2 12 ]-Vha'9 1 9 1 [■■■] = 

= m\ 2 \ha' [-..]. (28) 

The dots inside the square brackets stand for the terms that we have already 
calculated. This procedure can be generalized to give a series in h. Sadly, we 
cannot say much about the exact value of the integrals, so we write it this way 

oo 

If = E h 4 ( N2 -^ +k a^ N2 -V +k D( 8N2 - s+2k \y) (29) 

k=0 

where 

(y) = det 1 ^ (iph-'a'-'R^YZF^) (30) 

and each is homogeneous of degree / in the coordinates and in general, in 
the elements of the matrix M. With this, we have exactly calculated the classical 
term and obtained the dependence on the physical magnitudes of the quantum 
series. 

On the other hand, this is only valid for TV > 1. If N — 1, no fermionic 
operator appears in the Hamiltonian and the integral yields null. This is to be 
interpreted as a consequence of the discrete nature of the spinnorial degrees of 
freedom. In other words, the 'volume' occupied by the Grassmann variables has 
null measure according to Berezin's rules of integration. The appropriate count 
of the fermionic degrees of freedom in this case is the sum over spin polarizations 
that we made for the free configurations (classical geometry). When the Hamil- 
tonian depends on the spin, then the integral does contain all the information. 
The classical configurations are included although they do not contribute. Both 
in the fermionic and the bosonic case, they lie on spaces that have zero measure 
if compared to the global phase space. 

In order to perform the bosonic integral, we expand the exponential sepa- 
rating the terms in the Hamiltonian with coordinate or momentum dependence. 
Symbolically we write 
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-H(p)-H(y) 



n=0 



-1^ 



H{p) + H{y) 



(31) 



Where p and y stand for any momentum or coordinate operator. We shall keep 
this notation from here on. Once the binomial is expanded, we are left with a 
series in the coordinates and the momenta with no particular ordering. Taking 
advantage of the canonical commutation rules we can expand each term as a 
finite series in h. The coefficients of the series will depend on the particular order 
of each term and on whether the momentum and coordinate operators do have 
the same indices or not. An example is 



p 2l f m 



disordered 



ordered 

+h 2 x c2 x p 



+ h x cl x p 

y 



21-1 Mm-1 



ordered 



+ 



ordered 



+ O (h 3 ) 



(32) 



with cl and c2 some constants. 

The integration is carried out inserting the expansion of the identity between 
the coordinate and momentum operators so as to transform them into numbers. 
This insertion need just be done once in each term thanks to the ordering. This 
way we make only one integration over the phase space and we shall be able to 
extract common properties of all terms. One of them would be integrated to give 



I b = J ^^^- 8 ^+ 2s+8 ^ +fc D(^ 2+2fc )(y)a s (ap 2 )'- s 6 s / 2 (% 4 ) m - s / 2 . (33) 
Where 

a = and b = — . (34) 

2h 167r 2 ftV 3 V ; 

The number of commutations that have been needed to order the operators is 

2s. We only consider even terms because the others are related to odd integrals 

that yield null. We shall ignore again the contribution of the identity, as a 

generator of U(N). It can be easily calculated later. Now we change variables 

p — > a~ 1 / 2 p and y — > b~ 1 / 4 y and take advantage of the homogeneity properties of 

the D-functions to simply get 

J dp'dy'D^ N2+2k \y')p' 2l - 2s y' im - 2s . (35) 

It is manifest from here the remarkable scaling property that the dependence on 
the physical magnitudes f3 and R is the same for all terms in the series indexed by / 
and m. This is quite fortunate since all the bosonic integrals are, by themselves, 
polinomically divergent. We had to make this expansion in order to get the 
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quantum series that disordered is hidden in the expansion of the exponential in 
Eq. fl3i~D. Now we can sum up the series over / and m before integrating. These 
are the integrals that should be finite, at least for physical reasons, because they 
are the coefficients of the quantum series. 

Namely, what we wanted was to expand the original exponential of the Hamil- 
tonian in terms of the h that appears in the commutators. We know that this is 
feasible because it is physically meaningful. Then we expanded the exponential 
and the binomials that appeared, and after that we expanded again each term 
as the series in % that we were looking for. As we have found that the physical 
behaviour of the terms depend only on the index s of the quantum series, we are 
allowed to sum back again every term with common s, that is, to undo the first 
two expansions. Therefore the partition function for a fixed N is 

ziN > 1} = g g g{Nj s k) tPR\ lS+lk e -Hiv (36) 

s=0fc=0 V a ) 

As we mentioned, this is not valid for N = 1. 

The way we have made the calculation seems to appear a little clearer if we 
separately compute the 'classical' s = 0, k = term. That amounts to ignoring 
the commutation rules from the beginning or equivalently truncating the series 
in Eq. fl32|) to the first term, i.e. the 0(1) in h. That way, we can reconstruct 
the exponential and get 

Z^rrL = J dp < ple-^det^M^e-^b > = 
J dpdp'dxdx' < p\e~ pH{p) \p' >< p'\x >< x\det l/2 M(x)e~ pH(x) \x' >< x'\p > = 

J dpe~^ v) J dxdet 1/2 M(x)e-^ {x) .(37) 

To finally arrive at 

oo -B— ( \ 

*222 = E e -^r I (rW) ^in^F^e-^ . (38) 

N=2 ly ■ J \k,X J 

We explicitly show all the indices for this simple case and again, we have ignored 
the N — 1 term. The iV! factor comes again as the remnant of the gauge symme- 
try along the minima of the potential. It is precisely the same result as in Eq. (^) 
and we can see that the integral that appears as a coefficient is finite and perfectly 
defined. Apart from the dependence on the temperature, the expression differs 
from the calculation with the free configurations in the absence of the volume in 
the final result. This is due to the fact that volume is a physical quantity that 
only has a complete sense for those configurations. Here, it is taken into account 
in the flat, divergent directions of the integral. There, we would have different 
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powers of the volume depending on whether we integrate over the DO-brane wave 
functions or their bound states. So even when we do not see powers of the nine 
dimensional volume, our computation is a partition function and as such it is 
not a linear function of Vg. The fact that these integrals include by themselves 
the effects of bound states can be seen by just inserting a (— 1) F operator. This 
would turn this function into the Witten index, which is precisely used to count 
such states. These states also include different dependences on the volume even 
before perfoming the sum in N, because they act as if certain directions along 
the classical phase space were frozen. 

Let us now complete the calculation by adding the N — 1 term, which exactly 
equals the integral over the fields related to the identity -the center of mass of 
the system. This is special in several ways. Not only is it independent of the 
fermionic variables, but neither does it depend on the bosonic coordinates. Since 
it is just a function of the nine bosonic momenta and the Gauss constraints are 
trivial, it does not receive any quantum correction. This is clear when we relate 
this special matrix configuration with a single object in eleven dimensions: one 
supergraviton on its light-cone. One single object is always free and its statistics 
is irrelevant. It was already computed in the previous section, so we simply write 



Z = 256 V 9 (2nRph) 



-9/2 



e -S + EE/W9)(f 

N=2 q=0 V ' ' 




(39) 



We have been able to gather the two quantum series -the one of the fermions 
and that of the bosons- into a single one because they are expansions in exactly 
the same parameter. One of the most amazing characteristics of the partition 
function here calculated is its peculiar dependence on h. It only appears in the 
center of mass contribution and in the exponential coming from the longitudinal 
sum. For all the N > 1 terms we have supposedly made a classical approximation 
followed by a quantum series that would correct it; nevertheless the fc's that we 
include in the commutators are exactly canceled by those that appear in the 
Hamiltonian. So we arrive at the bizarre conclusion that the 'quantum' series 
is, in fact, not quantum at all. Moreover, this series is an expansion in the 
same parameter as the series in the longitudinal momentum N, therefore, we 
are tempted to conclude that the true natures of both series are not clearly 
distinguishable. The function is not very different from its % — > limit, the only 
role of this constant seems to be to serve as a kind of chemical potential that 
measures the cost of adding longitudinal momentum to the system. 

This poses the important question of how to relate this result with the ap- 
proximations we have made in previous sections and eventually, how to take the 
limits that would lead us to the different string theories. Different classical limits 
will not appear in this theory just as it is usual in traditional quantum theories. 
According to them, specific quantum statistical properties were effects related to 
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the nature of the fundamental fields. However, in Matrix Theory we have seen in 
a previous section that the only way of including quantum statistics is to assume 
that they appear as a consequence of the residual gauge interactions. This con- 
nects with the surprising behaviour of the quantum series. That is why we shall 
see that all these effects come out because of the appearance of new degrees of 
freedom related to the off-diagonal terms. The classical limits are reached when 
those degrees of freedom decouple getting infinitely massive. 

Nevertheless, the decoupling of degrees of freedom is always obscure when 
we look at the partition function because physical magnitudes are related to its 
logarithm and derivatives thereof. Different contributions to, for example, the 
Helmholtz free energy are not summed up but multiplied inside the partition 
function. Therefore, if some of the effects are taken to be negligible, the whole 
partition function may tend to zero. The best way to circumvent this difficulty 
is, then, precisely, to calculate the Helmholtz free energy. It is, basically, the 
logarithm of Z, which carries the advantage of turning products into sums, and 
so, clarifying the discrimination of the contributions. 

Let us, for example, take the classical limit 

7 7 7 

^classical^other 

(3F = - In Z classical - In Z othcI (40) 

We know that the classical partition function diverges when h — > so that the 
other term must tend to zero in order for the complete function to remain finite. 
Looking at the expression for the free energy one can see that the non-classical 
term acquires an infinite energy while the classical one gets less and less energetic 
as we get nearer to the limit. 

It is clear that the factorization of the partition function is quite arbitrary 
as we can separate any two parts by putting the appropriate parameter and 
then decouple one of them by taking a certain limit for the parameter. This 
is a consequence of the fact that we can organize the degrees of freedom of the 
theory in many ways -with or without physical meaning- and that, indeed, the 
information that the partition function holds is global so that the behaviour 
of particular configurations may be, as we have seen, quite hidden. The only 
physical way of knowing which is the correct expansion in each case is to go back 
to the Hamiltonian, then decide which degrees of freedom we are interested in 
according to the limit we are choosing (classical, large distances, ...), and only 
after that should we integrate exclusively the appropriate configurations. Some 
limits will be calculated in next section. 

Let us not forget that the series is itself a limit in the sense that it is only 
a good expansion for (f3R) — > 0. Its validity is basically determined by the 
coefficients f(N,q). What we know is that the f(N,0) carry a 1/N\ factor with 
more N dependence coming from the Gaussian-like integrals and that the other 
coefficients come from a quantum series that should be well behaved. If we fix 
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the radius of the eleventh dimension, and look at the series at high temperature, 
we can see that the thermodynamical behaviour is completely determined by the 
coefficients f(N, 0) but it is clear that the quantum corrections are the only ones 
that are affected by variations of the temperature. The question of what is the 
mechanism through which the theory remains finite without any ultraviolet cut- 
off is unresolved yet, because it is related to the growth of the coefficients with 
N. This would tell us how the theory distributes its degrees of freedom. An 
analogous high temperature expansion can be obtained for perturbative strings 
by making a high temperature expansion for the free energy of each field in the 
string. This always gives a Laurent series for (3 — > with leading term f3~ d with 
d — 1 the number of open spatial dimensions for the Helmholtz free energy. The 
point is that the number of degrees of freedom per mass level grows too fast so 
as to get the high temperature series a non convergent one. Actually, no high 
temperature divergence appears because the canonical equilibrium gets broken 
at the Hagedorn temperature. This is the way modular invariance works as an 
ultraviolet cut-off at finite temperature in fundamental strings. 

On the other hand, for a fixed temperature, the series is not able to describe 
systems with large light cone radius. It seems to be adapted to be a good descrip- 
tion of type IIA String Theory. Indeed, one just has to make the substitution 
R — > g s \Ja'% to recognize the series in iV as a String Theory non-perturbative 
expansion with the exponential of the inverse of the string coupling as expansion 
parameter. Each term is further corrected by a weak coupling expansion, which 
represents the quantization around the different solitonic sectors. This quantum 
series recovers its natural parameter K. So this would be the partition function 
of type IIA String Theory including the non-perturbative effects. Nevertheless, 
there is no exclusively perturbative term in this series so that strings themselves 
do not seem to be accounted for. This is a consequence of the particular choice 
of reference frame that we have made in the eleven-dimensional theory which 
we began with. The length R we have been writing is not exactly the R n that 
is what is directly related to the string coupling but the same radius measured 
by a strongly boosted observer. The same happens to the temperature, it is an 
energy and so it is affected by Lorentz transformations. Anyway, this is not much 
worrying since we have calculated a scalar dimensionless quantity; that is, the 
partition function is the same no matter if it is calculated in the light-cone or 
in any other reference frame. This way, we can invert the boost and make the 
substitution 




(41) 



with 



R u 



(42) 



R 



oo 
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It may be useful to remind that when we say 'at rest' we mean that the observer 
measures the invariant radius of the eleventh dimension, the smallest, so that, 
somehow, he is at rest with respect to the compact circle. In this reference frame, 
the parameter is much bigger so that the series is useless except for extremely high 
temperatures. Therefore, in spite of being Lorentz invariant, the series is only 
adapted to describe light-cone objects. Now the relation to the type IIA string 
is more direct but still, we do not seem to have any string! The reason for this is 
that our light-cone calculation has 'integrated out' the string degrees of freedom 
and spreaded them over the DO-brane ones. The degeneracies that are counted 
by the series in the index N do not either come from strings nor DO-branes but 
rather from what in Solid State Physics would be called quasi-DO-branes . These 
are the original branes, but dressed by the closed type IIA strings, other D-branes 
and anti- D-branes and interacting through them. Those are the solitons to which 
the series refers and this is the reason why we do not have any zero mode. The 
conclusion is that our calculation corresponds to the partition function of type 
IIA strings, but it is adapted to a point of view different from the usual one. 



5 The field-theoretical approach 

In the previous section we have computed the partition function for the Matrix 
Model using the commutation properties of the matrices that describe the DO- 
brane dynamics. This section is devoted to a calculation of Z(j3) using the per- 
turbative expansion for the U(N) Super Yang-Mills Quantum Mechanics. This 
type of calculations are closely related to those done in || [TT| , where the one loop 
effective potential between two DO-branes is obtained. We will also compute the 
thermal degrees of freedom for the case of the Matrix Model on T 2 , using the 
dual description in terms of the D-strings dynamics. Comparing both analysis 
we will check T-Duality. 

It is assumed that, when studying the interaction of non-BPS states, the 
correspondence between Matrix Model and Supergravity calculations may be 
realized at finite temperature ]Tj§, [15|, This picture comes from the idea that 



the non-extremality of the D-branes would be included into an entropy, and it is 
closely related to the Supersymmetry breaking produced by the different statistics 
of fermionic and bosonic fields^]. 

We will work with the Euclidean version of the SYM action, that is (in the 
temporal gauge, A = 0) 



S 1 



dr Tr 



R 



OL 



167r 2 a /3 



+ 



R 



2na 



(43) 



3 In a recent paper closely related to this part of our work, A. Tseytlin has analyzed 
these ideas. 



18 



where the Euclidean time r is taken to be a circle of length (3 and we have set 
h = 1 again. The Y l are nine scalars fields and the are Grassmann fields of the 
transverse SO (9) rotation group. Both are matrices of the adjoint of U(N). 

We will compute the one-loop correction to the free partition function for the 
fields that appear in the Matrix Lagrangian, and compare to the results obtained 
in the previous sections. The physical meaning of this type of calculations will be 
explained in detail below. We will show how this description, like the calculation 
presented in [P7|1, only takes into account the thermal properties of the fields 
attached to the DO-branes without thermalizing the D-particle itself. 

As usual in quantum field theory, we can obtain the corrections to the partition 
function by using the Feynman diagram techniques. Expanding the interaction 
term of the action in the Path Integral functional we have 



Z{P,N) = Z ((3) 



1 + E 



■l) n 



n=l 
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Z n {(3) 



(44) 



where Z n stands for the n-vertex correlation functions. As an example, in a 
simple model with two Yukawa-coupled real fields we would have 



Zn = Jf[ d Xi d y:i dz k J I?0^(A0(x i )^(%)0(^O) n e-( S "- W+S "- [V ' 1 ) (45) 



and 



Z {P) = [det (K(x,f3,<P))}^~ [det (K(x, (3,^))}^ 



(46) 



with F the fermionic number. The above expressions may be trivially generalized 
to the momentum space. 

We will start by computing and studying the tree level contribution to the 
Helmholtz free energy by using the Schwinger representation of the field propa- 
gator. Following this method we write 



-Tr 



& ~{k 2 +m 2 ))t 



t 



(47) 



In our case we have massless fields in d — + 1 dimensions, then the trace in 
the above expression has to be taken over the Matsubara frequencies and over 
internal degrees of freedom which are the representations of the transverse 5*0(9) 
and the U(N) gauge group. As the result of a blind calculation, we would propose 



7 N 2 r°° -, 
-/3F 0/ 3=-— / dtr 1 
2 Jo 



90 3 0,— - 80 2 0, 1 



(3 2 J 



(3 2 J 



(48) 



Here the N 2 factor comes from the U(N) degrees of freedom, and the multi- 
plicative factors of the thermal modular functions reflect the number of bosonic 
(9) and fermionic (8) fields. Some comments about the previous expression are 
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needed. If we were studying a second quantized theory of objects living in eight 
dimensions we would expect the Helmholtz free energy to grow proportionally to 
the eight-dimensional volume. As we see Fq{/3) in (f!8f ) is volume independent. 
This fact shows that what we are really doing is to thermalize the internal degrees 
of freedom of the DO-branes, that is, the strings attatched to them. 

This type of situation has already been studied from the string point of view. 
However it seems that there are some subtleties to be taken into account to 
establish the possible connections. In [17] the thermal free energy of an open 



superstring gas in presence of a Dp-brane has been computed, obtaining for the 
special p = case 



N 2 r°° dt 



6* 



2nit 



MO. 



2mt\ 



/(*) (49) 



where f(t) is the partition function of the Type I string theory. If one takes the 
zero temperature limit one recovers the vanishing vacuum energy of a supersym- 
metric theory. This property have been used by Polchinsky [18j to show the BPS 
nature of the D-branes. If we want to compare this expression with fl48|) we must 
consider only the massless field content of the string spectrum giving 



8 3 0, 
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(3F m =°((3) = -— dtr 1 
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(50) 



Where the 8 corresponds to the number of degrees of freedom of the massless field. 
Remember that at this level the string spectrum is composed by a vector and a 
Majorana-Weyl fermion. Both belong in their respective SO (8) representation. 
However in ([481) we would have obtained a different result that could be written 



„ /v 2 r°° ( 

- (5FM = -(3F m =\(3) -—j Q dtr% 0, 



2nit 



(51) 



The last term in the previous equation prevents F (j3) to vanish at T = because 
of the 03(0, s) zero mode 0. This property of the Helmholtz free energy we have 
obtained for the Matrix Lagrangian seems to break the supersymmetric nature 
of the theory in such a way that the D-particles would become non-BPS states, 
but there are some details related to the gauge we have chosen that will restore 
our standard knowledge about D-branes. 

In the string calculation, the Dp-brane configuration is obtained by compact- 
ifying Type I string theory on a torus with radii Rp+x, ...,Rd-i which then are 
taken to zero. The free energy is then calculated integrating the momentum of 

4 If we neglect the zero mode of the fields, we see that both expressions vanish in the zero 
temperature limit, but they will do it in a different way. 
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the string coordinates with Neumann boundary conditions. When we do it we 
fix the D-brane at a given position in the Dirichlet direction; only allowing fluc- 
tuations on the D-brane world- volume. In other words this type of computation 
only takes into account the string's knowledge of the full d- dimensional space- 
time, that gives the eight bosonic and fermionic contributions in (f£9|). In the 
particular case of a D-particle we assume it to be fixed at a given point in space, 
and we compute the string contributions to F(/3). From the point of view of the 
SYM Quantum Mechanics that describes the D-particle dynamics, and then from 
the Matrix Model one, this kind of configuration must correspond to an election 
of a fixed background. However this is not enough. If we let all the D-particle 
coordinates fluctuate we would once more have nine bosonic contributions with- 
out connection to the string calculation. What we are forced to do is to assume 
that there is one non-fluctuating direction. To choose the frozen direction we can 
take advantage of the worldsheet conformal invariance of string theory that for- 
bids the vibrations along the longitudinal direction of the string. In our case this 
coordinate coincides with the straight line that connects the D-particles [JXTJ, O . 
Consequently, if we want to relate the Matrix Model computation and its string 
origin we have to freeze this direction. 

We can analyze this problem from a more technical perspective. As we said 
the theory that describe the D-particle dynamics is a SYM quantum mechanics 
coming from a dimensional reduction of the theory ind = 9 + ltod = + l. 
We can start our analysis taking the Arnowitt-Fickler gauge in the initial SYM 
quantum field theory, that is 

n,A\ = (52) 

where n M are components of a unitary vector. Here the gauge fields are functions 
of the full ten-dimensional space, then after the dimensional reduction we have 

n<,^ + 71^ = j = l,...,9. (53) 

Remember that the condition in (|52|) does not totally reduce the degrees of free- 
dom to the physical ones. So we must take another gauge condition to complete 
the reduction. Here we can choose the temporal gauge (A$ = 0), that reduces 
the relation in (^3j) to a constraint between the scalars fields Yj. This picture is 
what underlies the physical idea of freezing the vibrational modes of the string 
along the longitudinal direction, and corresponds to the picture described in the 
previous section for the hamiltonian formalism. 

We could have done the same analysis by choosing another gauge and intro- 
ducing the adequate ghost contribution. As in || it is possible to choose the 
background field gauge, in the special case of the coinciding D-brane configura- 
tion. In this case the spectrum of the theory is composed of eight on-shell bosons, 
after considering the ghost contribution, eight fermion fields, and the interactions 
between them are the same we have in the temporal gauge. 
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To summarize, we have then shown that, for the d = + 1 gauge theory 
describing the dynamics of our system, the temporal gauge does not completely 
reduce the degrees of freedom to the physical ones. Another reduction is then 
needed. This constraint has to be taken from the dimensionally-reduced theory 
we started from. In the case of the D-particle dynamics it is its string vibration 
origin what determines the reduction. 

Let us come back to the volume dependence of (|50|). From an exact M- 
theoretic point of view we should compute the contribution coming from the 
nine directions. In fact we would expect a linear volume dependence of the 
Helmholtz free energy. These contributions could be studied by assuming that we 
are separating the physical degrees of freedom by means of a Born-Oppenheimer 
approximation [O]. We are assuming that the phase space of the theory may be 
factorized in terms of the dynamics of the D-particles as objects moving in the 
nine transverse directions and their vibrational degrees of freedom. The eight 
transverse string modes are decoupled from the motion of the DO-brane which 
must be studied as we did in the previous section. Finally we may describe the 
full dynamics in terms of a phase space coming from the motion along the nine 
dimensions and characterized by the presence of the string vibrational degrees 
of freedom for each point in that space. This idea have been also proposed in 

Em. 
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relation with p-brane and black-hole physics 

Let us start a more detailed study assuming that the N DO-branes are located 
in a background characterized by 
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(54) 
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This election corresponds to a fluctuation of the D-particles orthogonal to the 
Y\ = Xi fixed plane. This configuration completely decouples the Y\ field from 
the action. On the other hand we consider the transverse vibration of the massless 
fields connecting the branes. We should do that taking the following parameter- 
ization 



{Y k )\ = x k S) + A$ 



(55) 



with 0* G U(N). This corresponds to coinciding DO-brane positions, preserving 
the full U(N) gauge symmetry. This configuration exactly coincides with the 
string approach to the thermal strings stretched between N D-branes in W7 . 

It is possible to generalize our study to the separated branes configurations. 
As it is well known, separating the D-branes breaks the U(N) symmetry down 
to U(Ni) x [/(AT 2 ). It is easy to see how this breaking occurs. As an example we 
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can look at the N = 2 case. Taking the branes in different positions in the frozen 



direction, we have 



Y-i 



X\ 

2 2 



and Yi 
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(56) 



where T a are the group generators. One can easily obtain the mass correction to 
the scalar fields, coming from the interaction term in the action which is 



Tr [Yi, Y^' 



a 1 {X 1 - X 2 ) 



+ 



[2/2 



(57) 



that is diagonalized in terms of the massive fields <pj and </>?. This property was 
firstly pointed out by E. Witten in f20fl. Finally the free energy reads 



(3F(f3) 



(5? 



The same arguments hold for an arbitrary iV case. If we separate this system 
into two objects with Ni and N 2 , the free energy will have the form 



N 2 N 2 



(59) 



Where G\(0) and G 2 (/3) are the functions of the temperature corresponding to the 
thermal degrees of freedom of the fields living attached to iVj block of D-particles. 
Following this procedure it is possible to change the N 2 dependence in ( |50"D into 
a linear one on N, that corresponds to a U(1) N gauge group. The description we 
have shown in terms of the string physics included in the DO-brane dynamics may 
be done without any reference to it. As we have said for the overlapping D-brane 
configuration, here we may also fix any other gauge condition and obtain the same 
results as with the temporal (Arnowitt-Fickler) gauge. The spectrum obtained 
pulling apart the particle positions is, as expected, gauge independent. It is worth 
noticing that, as in quantum field and string theories |21j, the interaction can be 
partially absorbed as a correction for the masses of the fields. 

Now let us carry out the computation of the one-loop correction to the free 
energy. As we did at tree level we have to decide how many string directions 
we take into account. Assuming the philosophy of the Born-Oppenheimer ap- 
proximation we have explained before, we will fix the corresponding longitudinal 
direction. Following |8], || we rewrite the action ( fffi) in the usual SYM form, 
which reads 



S E = X f drTr 



g 2 Jo 



yiyi + qTq + _ 



Y\ Y j 



Y\e 



with 



R 3 / 2 
2vra' 3 / 2 



(60) 



(61) 



23 



An adequate rescaling of the bosonic and fermionic fields with the SYM coupling 
g allows us to study the Feynman rules of this theory. It is also useful to explicitly 
write down the matrix fields in terms of the U(N) generators, obtaining^ 

S E = J*dr Y:Y: + 6% + 9 ^YMY}YjF abc F^ + g6^Y:6 d F c 



^cdb 



(62) 



To be complete we can write a formal expression for the n-vertex correction to 
the Helmholtz free energy, which is 



FG9) = ~logZoC9)-^log 
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(63) 



We show how to compute the functions Z n for a very simple case in fl45). In the 
case of the Matrix Model we shall write Z n {0) as 
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where I(Y*) and I(9 a , Y*) are the interaction terms appearing in the Matrix 
Lagrangian. 

Now we can explicitly write the Feynman rules of the Lagrangian in ( |62|) 
Scalar propagator = — wit hp = Po = —5— {n G Z) 



Fermion propagator = - 

and the interaction vertices are 
,2 
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(65) 



in both cases the total momentum of the incoming fields is zero. It is easy to check 
that the correlation functionals in (|4l)D vanish for an odd number of fermionic 
fields, leaving the perturbative expansion of Z(/3) to be a power series in g 2 . 
Finally we must take care of the discrete nature of the momenta when we compute 
the loop integrals. The one-loop correction we want to obtain corresponds to the 
diagrams in Fig. [l|. 

The purely scalar contribution is 

Zl calars {f3) = -336 g 2 F ade F ade 7 ^Y, 
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(66) 



5 Here we take the normalization Tr (T a T b ) = S ab 
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(a) (b) 



Figure 1: The one- loop Feynman diagrams involving scalar fields (a) and the coupling 
of fermions and scalar fields (b) 



The factor 336 appears because of the 8 bosonic oscillating directions, i.e. after 
freezing the longitudinal one. 

The contribution coming from the diagram involving fermionic fields gives 
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(67) 



where we have taken the trace over the fermionic degrees of freedom included 
into the jj. We are now able to study some properties of the one- loop correction 
to the partition function, which takes the form 



Z 1 ((3) = -Zr lars (P) + -Zt m \f3). 
Then the corresponding free energy looks like 

F{p) = ~togZ ()9)-itog[l + Z 1 ( ) 9)] = 



F (/3)--log 



i + g 2 



(N 2 -l)(N 2 -2) 



7T 3 



G 



(68) 



(69) 



25 



where G is a constant coming from (p5|) . The value of this constant can be 
obtained by a generalized-^ function regularization. In the zero-temperature case 
this contribution vanishes by dimensional regularization because of the fact that 
there is no dimensionful parameter involved. This means that we are dealing with 
a high temperature expansion whose relationship to the T = case is subtle. 

We can now check what happens if we consider the Matrix Model compactified 
on a torus. The dynamics of N DO-branes on T n is described by a system of N 
Dn-branes on the dual torus ]T|, f|, [EJ. In our case we start with the Matrix 
Model on S l (R) x (S l (L) x R 8 ), where the first S 1 stands for the compactified 
eleventh dimension. In this case the dual description is given by a system of D- 
strings wrapped around a cylinder of length S = 1/L. The theory that describes 
this system is a d = 1 + 1 SYM quantum field theory coming from dimensional 
reduction from d = 9+1. Because of the appearance of the £ = 1/L multiplicative 
factor in the dual action El we can redefine the coupling constant of the theory 



in such a way we can write 

glai = f = g 2 x (70) 

where g is the effective coupling constant in the Kaluza-Klein sense. 
In this case the gauge condition in fl52|) reduces to the following form 



n A° + mA\ + nj Yf = (71) 

then the problem of gauge fixing in this case is exactly the same as the one we 
have analyzed before. After taking the temporal gauge we have to choose the 
additional condition coming from the residual constraint that relates the eight 
scalar fields Yj and the A\ gauge field. Initially we have two possible options, 
either we can set the gauge field to zero, assuming that there are no propagating 
modes on the D-string, or we can stop the string vibration in one direction. In 
the latter case we will have a purely Neumann string mode propagating on the 
D-string and interacting with seven scalars. In both cases we will have eight 
bosonic and fermionic degrees of freedom. At the tree level the free energy is the 
same as in (|50|), except for the Kaluza-Klein modes contribution 



-/3F s (/3) = -— / dtr 1 
2 Jo 



8 , 3 0,M - 8 , 2 o, 2 " ( 



fc(o,^-) 



P 2 J \ 1 p 2 

8N 2 f°° , * ( 2mt\ n ( 2mt\ 
= ~ — L ^(0.^3(0,^). (72) 

As we can easily see, we recover the free decompactified limit, that is the Matrix 
Model result, when we take the dual radius £ to zero. 

It is now easy to compute the one-loop corrections to the Helmholtz free 
energy for the compactified case. The interaction terms we have are independent 
of the gauge reduction and exactly map into the open Matrix-Model ones. We 
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obtain 



a 2 B 3 

Z scalars^ = _ 33 g pade pade^ £ ^ 

s (2vr) „, m , 



-336 



(iV 2 - l)(iV 2 -2) 



EEt 



n 2 +(m§) / 2 + (s|)' 
1 1 



n 2 + (m|) Z 2 + (sf ) 



• (73) 



To check T-duality in this expression, we have to rewrite the dimensionless cou- 
pling g in terms of the initial Matrix parameter, the coupling constant g and 
finally take the limit £ — > 0, which reduces this to (p6|) . The same arguments 
hold for the fermionic contribution (0). We can give a general argument for the 
n- vertex function in fl64}) . In this case, taking the limit of the open Matrix Model 
(that corresponds to the dimensionally reduced D-string system), one gets 



Z n 



oc 



(2vr)- 2 " ^ 




(74) 



This parameter is explicitly independent of the compactification radius, therefore 
it is self-dual and then coincident with that of the already open Matrix Model. 

The expression in fl64|) corresponds to the term with g 2n (or g 2n , when we 
compactify one dimension) in the perturbative expansion of the free energy 



F s (/5)=F s (/5)--^log 



oo ~2na3n 

1 + V^— G 



n=0 




(75) 



where G n is the function that comes from the diagram of order g 2n . The regu- 
larization arguments given for (|69D hold here too. This function depends on the 
Kaluza-Klein modes of the compactified dimension but, in the small £ limit, it 
only takes into account the corresponding zero mode. The behaviour of the free 
energy of a D-string system ( Eq. (|75"D) allows us to explicitly check T-duality 
order by order in perturbation theory. In fact if we computed a general n-loop 
contribution we would see that its dependence of the coupling constant, the tem- 
perature, and the length of the compact dimension would map, after taking the 
corresponding limit, into the d = + 1 SYM computation of the same contribu- 
tion. This property is not surprising because of the string origin of the T-duality 
of the D-brane dynamics. 

We can now go back to the physics involved in the Born-Oppenheimer ap- 
proximation. The complete partition function of the finite-N Matrix model has 
to be obtained multiplying the internal (stringy) degrees of freedom and those of 
the translational D-particle dynamics. In terms of phase spaces we can express 
the complete phase space as the direct product of the translational and internal 
ones. This assumption amounts to the following form for the partition function 



3orn-Oppcnhcimcr 



(P) = Z^iP) x Z IntcI ™\P). 



(76) 
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Finally, let us mention that adding the internal partition function to one apropiate 
term in the sum in (|17|) we recover the desired TV 2 dependence coming from the 
internal unbroken U(N) SYM degrees of freedom which are necessary to fit the 
entropy of a given iV-parton in (|17|) as that of a Schwarzschild black hole of the 
model in ||, and obtain the Bekenstein-Hawking area law. 

6 Conclusions 

Matrix theory shows several different aspects when one looks at its thermody- 
namical behaviour. When the distances between objects are large, they can be 
described by a semiclassical expansion around what can be later interpreted as 
the Kaluza-Klein modes of a supergraviton in eleven dimensions. This is, in fact, 
the low energy limit of M-theory and it is the only case in which we have been 
able to take the radius of the eleventh dimension to infinity and see what are 
the predictions of the model for a completely open universe. We have seen that, 
as maybe one could have expected from the beginning, it precisely corresponds 
to the results obtained for eleven- dimensional supergravity, defined in an open 
space. 

One interesting point would be the connection between our work and the fact 
that the semiclassical Born-Oppenheimer approximation for the supergraviton 
modes justifies their election in || as candidates to describe black holes. The 
contribution of their center of mass does not spoil the arguments exclusively 
based on the Super Yang Mills entropy. Another point of view is taken in 0, 
where the model supposes that a black hole consists of a gas of distinguishable 
D-particles with Boltzmann statistics. In this case, the separation between fast 
and slow modes of the system helps us justify the dependence of the entropy with 
the number of D-particles. Indeed, the field theory gauge symmetry is broken to 
U(1) N and it seems that S oc N may not be taken as an assumption. It keeps 
being unclear how the discernibility is acquired in this model. Maybe one could 
get some insight about this possible connection from the calculation in section 
four, where the whole gauge group is included. The question is that we do not 
really see how the black hole enters into the stage in our picture. 

In section four we have tried to calculate some global aspects of the parti- 
tion function including interactions and quantum effects. We have seen that the 
consequence of taking canonical Heisenberg commutators between the coordinate 
fields and their conjugate momenta is the appearance of a series whose natural 
parameter is not % as usual but the light-cone radius. In fact, the parameter is 
the same as that of the series coming from the sum over the possible values of p + . 
We see this as a confirmation of the idea that this theory is quantum in a broader 
sense that the presently known quantum theories and that String Theories and 
Supergravities come out as classical limits of a moduli space. 

It is known, both in field and string theory, that the corrections to the free 
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theory coming from interactions can always be seen as the appearance of effective 
masses that correct the original ones. We wanted to see that in our case, but 
to do it, we needed to make a Born-Oppenheimer approximation. It consists in 
separating the 'movements' of the system in fast and slow ones. That is what 
we made in section 5. We suppose that we can divide the system in clusters of 
supergravitons with weak interactions among them. Each of these objects has 
low energy so we can separate the fields that propagate along their world-line 
(fast movements) and the displacements of their centers of mass through the 
target space (slow ones). The world-line calculation has been performed and 
represents a first attempt to include interactions as well as quantum effects that 
do not have statistical nature. We have obtained a high temperature expansion 
that corresponds to the series in the Super Yang-Mills coupling constant. This 
calculation and the equivalent one using T-duality and D-strings can be used to 
relate this approximation with one loop expansions in string theory like those 



made in [17]. They can also serve as a finite temperature check of T-duality, that 



holds to all orders in perturbation theory, as expected. 
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